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Abstract 

We compute the Isgur-Wise function using heavy quark effective 
theory formulated on the lattice. The non-relativistic kinetic energy 
term of the heavy quark is included to the action as well as terms 
remaining in the infinite quark mass limit. The classical velocity of 
the heavy quark is renormalized on the lattice and we determine the 
renormalized velocity non-perturbatively using the energy-momentum 
dispersion relation. The slope parameter of the Isgur-Wise function 
at zero recoil is obtained at /3 = 6.0 on a 24^ x 48 lattice for three 
values of mq. 



1 Introduction 



In the determination of the CKM (Cabibbo-Kobayashi-Maskawa) matrix el- 
ement \Vcb\ from experiment of the exclusive decay B D^*^^ the heavy 
quark symmetry plays an essential role, because the universal form factor 
^{y ■ v'), so called the Isgur-Wise function, is normalized to be one at zero 
recoil in the infinite heavy quark mass limit. However the differential decay 
width for this process disappear at the zero recoil point, so that one needs 
to extrapolate the experimental data to this point. The lattice computa- 
tion of the slope ^'(1) enables this extrapolation in a model independent 
way. On the lattice Bernard, Shen and Soni [|T| and the UKQCD Collabora- 
tion PI have computed the Isgur-Wise function using the propagating quarks 
(Wilson fermion and clover fermion respectively) for the heavy quark. Since 
the heavy quark mass should be below the lattice cutoff in this framework, 
present simulations are performed around the charm quark mass. An alter- 
native way of treating the heavy quark is to use the heavy quark effective 
theory (HQET) on the lattice. Mandula and Ogilvie formulated the Lattice 
HQET and applied it to a calculation of the Isgur-Wise function]^. Their 
simulation, however, suffers from much noise so that the extraction of the 
ground state seems to be difficult. This problem is well known in the calcula- 
tion of the heavy-light decay constant using the static approximation where 
one is forced to use some method to enhance the signal of the ground stated 
or to introduce O^I/uiq) terms in order to decrease the statistical noise 0. 

Our calculation is based on the lattice HQET keeping a part of the 
0{l/mQ) correction terms. The inclusion of the kinetic term of the heavy 
quark reduces the statistical noise significantly and enables us to extract the 
ground state reliably. We calculate the Isgur-Wise function for three values 
of the heavy quark mass and discuss ©(l/mg) effect on it. 

The plan of the paper is as follows. Section 2 contains the formulation of 
the lattice HQET including 0{l/mQ) corrections. In Section 3 we describe 
our numerical simulations and discuss the results. In this section, after brief 
account of simulation parameters, we first mention about the renormaliza- 
tion of classical velocity. Then the calculation of the Isgur-Wise function is 
described. Sec. 4 is devoted to the conclusion. 
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2 HQET on the lattice 



The action of the lattice HQET including ©(l/mg) corrections is 



X) 



-iv-D — 



2mQ 



h{x) 



(1) 



where v,. 



V, iv^) is a four- velocity of the moving heavy quark. 



is 

a lattice covariant derivative and F^i, is a chromo- magnetic (electric) field. 
This action is a generalization of the non-relativistic lattice QCD (NRQCD) 
for finite velocities. The first term corresponds to the action in the infinite 
mass limit used by Mandula and Ogilvie and the second term describes 
the O^l/mq) corrections. The heavy quark field Q{x) is expressed in terms 
of h{x) as 



^trngv-a 



h{x) 
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— ^cr F 



D). For simplicity we neglect the spin dependent 
from the action and the 0{l/mQ) correction 



where = 
interaction term 

term from the field redefinition, i.e. the second term in eq. (|]). It should be 
noted that the ©(l/mg) corrections are not included completely with this 
approximation, but the kinetic term describes the motion of the heavy quark 
inside the meson and gives a major ©(l/mg) effect. 

The heavy quark propagator is obtained by solving an evolution equation 
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where we used the equation of motion in the infinite mass limit 

{v^Di - w-D)h{x) = 

to remove the term from the evolution equation. Dih{x) and Dfh{x) are 
defined as 

Ui{x)h{x + i) — U-{x — i)h{x — i) 



DMx) 
DfKx) 



Ui{x)h{x + + Uj{x - i)h{x - i) - 2h{x). 



(6) 
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The cost for the computation of this deterministic equation is much less than 
one of obtaining the propagator for the Wilson fermion using some iterative 
solver. A parameter n is introduced to stabilize unphysical high frequency 
modes. Considering the evolution equation for free field, the stability condi- 
tion to assure the convergence of the equation is 

I 1 - -Ho{k) I < 1 (for any ), (7) 




where ki = sin ki and ki = 2sin y. According to the definition of derivatives 
(H), A; in the last term of (^) should be k. We, however, use k instead of k 
because the effect of this term is small for small v, and for convenience in later 
calculation. Similar condition is required in NRQCD [^], but in the present 
case it depends on the heavy quark velocity as well as rriQ and n. Because 

of the presence of term -^v ■ k high frequency modes of residual momentum 
do not converge with time evolution for a large value of velocity unless n is 
taken to be sufficiently large. On the other hand the 0{a) effect becomes 
larger for larger n, then much large n is not preferable. Fixing n for a value 
of mq, magnitude of velocities are restricted below a certain value. For the 
values of parameter used in the present simulation, J2i is restricted to be 
smaller than 0.25. 

Multiplying the projection operator to G{x), the heavy quark propagator 
in 4-spinor representation is obtained as 

Sgix) = G{x) ® (9) 

The Isgur-Wise function is extracted from the three point correlation 
function 

{Q\B{v')Vi{v' ,v)B{v)^\Q) (10) 

where Vi{v',v) is a temporal component of the vector current hyfj^h^ and 
B{v) is a local interpolating field for an initial or a final meson state qj^hy 
for which we take the same mass for both initial and final mesons. The three 
point correlation function is related to the Isgur-Wise function ^ {v ■ v') as 

Xs Xf 
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Xs Xf 

OC TUB ^{v ■ v'){v + v')o 

^^-Ef{tf-ts)-E^{ts-U) ^fQ^ > > t,) (11) 

where Xf = {xf,tf), Xg = {xg, tg), Xi = (0, ti) and the component proportional 
to (f — f')o is neglected since its form factor is ©(Aqcd/'^q) and {v — v')o 
is small itself in our velocity region. Taking the following ratio of the three- 
point functions the exponential factor and the renormalization factor for the 
vector current cancel and we obtain 

73 fj. J. J. \ Gy^y' (t f , tg, ti^Gy' ^y(t , tg, ti^ 

^V,v'[tf,ts,ti) = — — - — — T- - -1- 

We use this relation for the calculation of the Isgur-Wise function. 



3 Simulations and results 

We used 120 gauge configurations of 24^ x 48 lattice at /3 = 6.0 in the 
quenched approximation. Each configurations are separated by 2,000 pseudo- 
heat bath sweeps after 20,000 sweeps for thermalization. For the light quarks 
we used the Wilson fermion with hopping parameters 0.153 and 0.155, and 
obtained propagators under the periodic and Diriclet boundary conditions for 
spatial and temporal direction respectively. The critical hopping parameter 
and the inverse lattice spacing determined from rho meson mass are Kc = 
0.156986(50) and a"^ = 2.215(64) GeV. 

Heavy quark masses and stabilization parameters we used are 




where tuq = 1.8 nearly corresponds to the bottom quark mass. As mentioned 
above, heavy quark velocities should satisfy the stability condition (H) which 
leads J2i < 0.25 for our sets of mq and n. In Table 1 we tabulate our sets 
of V and v' used in the simulation. 
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The mean-field improvement is applied when we compute the heavy 
quark evolution equation. Link variables are altered as U^{x) U^{x)/uq 
where Uq = { \Upiaq. Y^'^- Upiaq. is the product of link variables along plaque- 
tte. We used the value Uq = 0.8776 measured on our configurations. 

3.1 Extraction of the renormalized velocity 

Before we proceed to the calculation of the Isgur-Wise function we mention 
about the determination of masses and velocities of the heavy-light mesons. 
The velocity of the moving heavy-light meson vr could differ from the one 
of the bare heavy quark v for the lattice HQET because of the violation of 
the space-time 0(4) symmetryj^, |^. We extract the 'renormalized' velocity 
from the simulation results nonperturbatively using the dispersion relation 
of the meson. 

The heavy-light meson obeys the dispersion relation 

p° = + ^ + i - i^-^'r) } (13) 

where p'^ is the meson momentum which is decomposed as p'^ = mpv^ + k'^. 
k' differs from our definition of the residual momentum k,, which is defined 
through — mqv^. Using a naive mass relation 

mp = mQ + Eq (14) 

where Eq is the binding energy, k'^ is written as k'^^ = k^ — E^v^. Substituting 
this relation in above dispersion relation, we obtain 

\ 2mp^v^' j V mpv^^j 2mpv° \ \^ ' J 

(15) 

From now on, we write E{k) instead of k^, and attach subscript 'R' to the 
classical velocity v for distinction from ones as the input parameter. 
On the lattice, above dispersion relation is rewritten as 

^-lEik)-Eik=0)] = 1 - 1% . (k-^^k) + (k-C4- kf] ] , 

(16) 
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where 

In Fig. |I]we plot the asymmetry of the energy — [e~'^(''')+'^(°)— 6^^*^^^-'+'^*^°'']/ 
from which we extract the 'renormalized' velocity vr using the relation 

= — H- ■ sm ki Tj2 sm — . (18) 

For Eq and mp in the right hand side of (plSf ), we use the values obtained 
at = 0. For the value of v^, bare one is used. These substitution should 
be justified by the independence of these values on the velocity. Obtained 
values of the 'renormalized' velocity are listed in Table ^ and slightly larger 
than 'bare' values v , but in most cases consistent with v in present error 
level. Ratios vr/v are almost independent on v as shown in Fig. ^. 

The heavy-light meson masses can also be obtained similarly using the 
relation 

E{k)+EiO) I E{~k)+E{0) 1 „7„ \ 

1-- — = ^—(k -{~k)A (19) 

which we plot in Fig. ^. The obtained values for mp are summarized in Table 
1^. In Table ^, we also list the results for binding energy of heavy-light meson, 
Eq, determined from correlation function at zero residual momentum using 
equation (|l^. In eq. (pT]), we again use the values at -u = for Eq and nip. 
Both of mp and Eq are almost independent of the velocity and well satisfy 
the naive mass relation mp = mg + Eq. 



3.2 The Isgur-Wise function 

In the calculation of the three point function we set = 8 and tg = 20 where 
it seems to reach the ground state of initial meson moving with velocity v. 
Fig. ^ shows the behavior of Rvy(tf,ts,ti) as functions of t/ for several sets 
of V and v' for mQ=1.8 and k=0.153. Signal is rather clean and we observe 
clear plateaus beyond t/ ~ 24. Then we take tf =24-27 as a fitting interval. 
For the largest mass mQ=5.0 signal becomes so noisy and we cannot identify 
a clear plateau. Nevertheless we treat them in the same manner as mq = 1.8 
and 2.5 with tf =24-25. 
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Our results for the form factor ^(f/j-f^) at k=0.153 are shown in Fig. ^ for 
mQ=1.8, 2.5 and 5.0. Error bars in the horizontal direction are coming from 
the statistical uncertainty in the determination of the renormalized velocity. 
Statistical errors increase with heavy quark mass as expected. In our v ■ v' 
region the form factors have almost linear behavior where we extract the 
slope of i{vR ■ v'pi) at vr ■ v'^=l, which is usually denoted as — p^, from a fit 
of our data to the the form 

i{vR-v'^) = 1 - p\vn-v'j, - 1). (20) 

for each rriQ and k are given in Table ^ with the values extrapolated to 
the critical hopping parameter for the light quark. mQ=1.8 roughly corre- 
sponds to the B meson mass where our data should be compared with the 
experimental data obtained by CLEO Collaboration 1 10 1. 



p^ = 0.70 ±0.17 this work (21) 
= 0.84 ±0.15 CLEO (22) 

Our data is consistent with that of CLEO considering the large statistical 
uncertainty. We also compare our results with other lattice results obtained 
using the propagating heavy quarks for which the heavy quark mass is around 
the charm quark mass 

p2 = 1.24 ± 0.26 (stat) ±0.33 (syst) [1] (23) 
= 0.9t2(stat)t^(syst) §. (24) 

These are again consistent with each other in their large statistical uncer- 
tainty. 

In order to see the ©(l/mg) effect we extrapolate our p^ to the vanishing 
l/mg limit as shown in Fig. ^ p^ becomes smaller when we approach the 
infinite mass limit. We need, of course, more statistics to quantify the size 
of the 0{l/mQ) effect and to extrapolate our data to the charm quark mass. 

4 Conclusion 

We calculated the Isgur-Wise function near the zero recoil limit using HQET 
including the kinetic term of the heavy quark. Spin-independent 0(1/ tuq) 
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effect is included in this approximation. The obtained value for the slope 
parameter is consistent with the experimental value, but the statistical 
error is still large. A large amount of this statistical error is due to the 
large statistical uncertainty in the determination of the renormalized classical 
velocity and it originates in extraction of binding energy of heavy-light meson 
for finite residual momenta. It is essential for reducing the statistical error 
to improve the signal for finite momenta using the smearing techinique for 
example. 

It is important to include the ©(l/mg) terms completely and to quatify 
the size of violations of the heavy quark symmetry in several semi-leptonic 
form factors. The lattice HQET is also applicable for studies of heavy to 
light semi-leptonic decay {B plv etc.) which can be used for extraction of 
Vub from the exclusive decay width in a model independent way. 
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TABLES 



V { initial ) 


v' [ nnal ) 


(0, 0, 0) 


( 0, 0, Vf ) 

Vf = 0.05,0.10,0.15,0.20,0.25 

( 0, Vf, Vf ) 

Vf =0.05,0.10 


( 0, 0, V, ) 

Vi = 0.05,0.10,0.15,0.20,0.25 


(0, 0, ±Vf),iO, Vf, 0) 

Vf = 0.0, 0.05, 0.10, 0.15, 0.20, 0.25 
( 0, Vf, Vf ),( 0, Vf, -Vf ) 

V/ = 0.05,0.10 


( 0, Vi, Vi ) 

Vi - 0.05,0.10 


(0, Vf, 0)(0, 0, -Vf) 

Vf = 0.0, 0.05, 0.10, 0.15, 0.20, 0.25 
(0, Vf, ±Vf ),{0, -Vf, -Vf ) 

V/ = 0.05,0.10 



Table 1: Initial and final heavy quark velocities used in our simulation. 
When we take the ratio Ryy{tf,ts,ti) rotationally equivalent combinations 
of V and v' are used. 
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K = 0.153 


bare velocity 


mq 


0.05 


0.10 


0.15 


0.20 


0.25 


1.8 
2.5 
5.0 


0.054(04) 
0.056(04) 
0.059(05) 


0.105(05) 
0.108(06) 
0.112(08) 


0.157(07) 
0.161(08) 
0.166(12) 


0.209(10) 
0.216(12) 
0.222(20) 


0.263(14) 
0.271(17) 
0.279(31) 


K = 0.155 


bare velocity 


mq 


0.05 


0.10 


0.15 


0.20 


0.25 


1.8 

2.5 
5.0 


0.052(06) 
0.054(06) 
0.055(07) 


0.102(07) 
0.104(07) 
0.107(10) 


0.152(09) 
0.156(10) 
0.158(15) 


0.204(12) 
0.208(14) 
0.210(23) 


0.255(16) 

0.262(21) 
0.261(35) 



Table 2: Results for renormalized velocities for each k and mq. These values 
are obtained from simulation data nonperturbatively using the dispersion 
relation. 



K = 0.153 


bare velocity 


mq 


0.00 


0.05 


0.10 


0.15 


0.20 


0.25 


1.8 


2.30(06) 


2.28(06) 


2.26(07) 


2.24(09) 


2.21(12) 


2.17(16) 


2.5 


2.98(08) 


2.94(09) 


2.90(11) 


2.85(15) 


2.78(20) 


2.68(28) 


5.0 


5.39(21) 


5.31(24) 


5.17(33) 


4.99(51) 


4.81(82) 


4.72(134) 


«: = 0.155 


bare velocity 


mq 


0.00 


0.05 


0.10 


0.15 


0.20 


0.25 


1.8 


2.30(07) 


2.29(08) 


2.27(09) 


2.26(12) 


2.24(15) 


2.20(21) 


2.5 


3.00(10) 


2.97(11) 


2.93(14) 


2.89(19) 


2.83(25) 


2.74(35) 


5.0 


5.45(25) 


5.40(29) 


5.28(40) 


5.16(61) 


5.09(100) 


5.20(169) 



Table 3: The kinetic mass of heavy-light meson for each mq and k obtained 
using the dispersion relation. Obtained values are almost independent of 
velocity. 
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K = 0.153 


bare velocity 




0.00 


0.05 


0.10 


0.15 


0.20 


0.25 


1.8 
2.5 
5.0 


0.531(04) 
0.533(04) 
0.533(06) 


0.530(04) 
0.532(04) 
0.534(06) 


0.529(04) 
0.532(05) 
0.535(08) 


0.528(04) 
0.531(05) 
0.537(10) 


0.527(05) 
0.530(06) 
0.539(13) 


0.525(06) 
0.528(08) 
0.544(19) 


«; = 0.155 


bare velocity 


mq 


0.00 


0.05 


0.10 


0.15 


0.20 


0.25 


1.8 
2.5 
5.0 


0.497(05) 
0.500(06) 
0.500(09) 


0.497(05) 
0.499(06) 
0.502(09) 


0.496(05) 
0.499(06) 
0.503(10) 


0.495(06) 
0.498(07) 
0.505(12) 


0.493(06) 
0.496(08) 
0.509(17) 


0.491(07) 
0.495(10) 
0.514(25) 



Table 4: The binding energy obtained from correlation functions at zero 
residual momentum for each velocity, mg and k. Obtained values are almost 
independent of velocity. 



K 




1.8 


2.5 


5.0 


oo 


0.153 


0.71(09) 


0.72(12) 


0.77(30) 


0.78(30) 


0.155 


0.70(13) 


0.71(18) 


0.77(47) 


0.75(46) 




0.70(17) 


0.69(24) 


0.77(66) 


0.73(67) 



Table 5: Results for = —^'(1) for each mg and k,. They were obtained 
by fitting 1 — p^{yR-v'^ — 1) with simulation data. Results at mq = oo are 
obtained from extrapolation in l/mg and the values at the critical hopping 
parameter ( Kc = 0.156986(50) ) are obtained by linear extrapolation in 
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Figure 1: Extraction of renormalized velocity at mq = 1.8 and k = 0.153. 
It shows -[e-^W+-^W-e-^(-'=)+^(°)]/2 vs sm{k,)-{Eo/mpv^^)2sm{kj2) 
for V = (0, 0, 0.10) and (0, 0, 0.20), where we use the obtained values at = 
for Eq and mp and bare values for v*^. Solid lines represent linear fits from 
which we extract the renormalized velocities. 



Figure 2: Sv/v for itlq = 1.8 and k — 0.153. 5v = vr — v where the values 
of vp are listed in Table 2. 5v/v are independent of v and slightly larger than 
zero, but consistent with zero in present error level. 



Figure 3: 1 - [e-^('=)+^(o) + e--E^(-fc)+^(o)]/2 as a function of (2 sin(A;^))2 for 
V — (0,0,0.20) at tuq = 1.8 and k = 0.153. Solid line represents linear fit 
from which we extract the kinetic mass of heavy-light meson. 



Figure 4: Ryy{tf,ts,ti) as functions of tf for several sets of v and v' at 
niQ — 1.8 and k, — 0.153. For all combinations, v's have the same absolute 
values as v and opposite spatial directions to v. 



Figure 5: The form factors at ttiq = 1.8 (a), 2.5 (b) and 5.0 (c) 

where n — 0.153 for all cases. Figures show the almost linear behavior of 
^{vr-v'^) in our velocity region. Slope of these functions, p^, are extracted 
by fitting to the form 1 — p^{yR-v'ji — 1). 



Figure 6: Extrapolation of to the vanishing l/mq limit k — 0.153. 
The values of p^ are listed in Table 5. 
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